Excitonic BCS-BEC crossover at finite temperature: Effects of repulsion and electron-hole mass 
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The BCS to Bose-Einstein condensation (BEC) crossover of electron-hole (e-h) pairs in optically excited 
semiconductors is studied using the two-band Hubbard model with both repulsive and attractive interactions. 
Applying the self-consistent ?-matrix approximation combined with a local approximation, we examine the 
properties of a normal phase and an excitonic instability. The transition temperature from the normal phase 
to an e-h pair condensed one is studied to clarify the crossover from an e-h BCS-like state to an excitonic 
Bose-Einstein condensation, which takes place on increasing the e-h attraction strength. To investigate effects 
of the repulsive interaction and the e-h mass difference, we calculate the transition temperature for various 
parameters of the interaction strengths, the e-h particle density, and the mass difference. While the transition 
temperature in the e-h BCS regime is sufficiently suppressed by the repulsive interaction, that of the excitonic 
BEC is largely insensitive to it. We also show quantitatively that in the whole regime the mass difference leads 
to large suppression of the transition temperature. 

PACS numbers: 71.10.Hf, 71.35.Lk, 03.75.Nt 



I. INTRODUCTION 

The Bose-Einstein condensation (BEC) of excitons in 
solids is one of the most interesting subjects in condensed 
matter physics. 1 In this phenomenon the attractive Coulomb 
interaction between electrons and holes plays an essential role 
in forming bosonic electron-hole (e-h) bound pairs. This in- 
teraction and the bosonic nature strongly affect the statistical 
and thermodynamic properties of e-h systems. 

In e-h systems realized in photoexcited semiconductors, 
therefore, various remarkable properties are expected to de- 
pend on e-h density, temperature, etc., and they have been 
investigated extensively both experimentally and theoreti- 
cally^ In particular, metal-insulator transitions have at- 
tracted interest for many years: One is the exciton Mott transi- 
tion from an exciton or biexciton insulating gas phase to an e- 
h plasma (normal) metallic phase with increasing e-h density. 
Another, which is our main target in this work, is the transition 
from the normal phase to an e-h pair condensed one. At low 
e-h density, strongly bound e-h pairs are expected to undergo 
BEC as an exciton gas at cryogenic temperatures. 4,5 On the 
other hand, at high e-h density where the mean interparticle 
distance is shorter than the exciton Bohr radius, weakly bound 
e-h pairs should behave like the Cooper pairs in conventional 
superconductors at sufficiently low temperatures, that is, the 
Bardeen-Cooper-Schrieffer (BCS) state of e-h pairs,-'? There 
have been many experimental attempts to observe such a co- 
herent condensation of e-h bound pairs in solids such as CU2O, 
CuCl and GaAs. 8 ' 9101112 At present, however, definitive evi- 
dence is still lacking in experiments. 

Since the e-h recombination is usually much slower than the 
intraband relaxation, we focus on a quasi-thermal-equilibrium 
state of the e-h system. In this situation, the crossover prob- 
lem 13 i 14 i 15 i 16 between the BEC of excitons and the e-h BCS- 
like state is fascinating, especially from the viewpoint of the 
difference from the BCS-BEC crossover in superconductors 
or trapped atomic Fermi gasesiiii&iS^ We believe that re- 
vealing such differences will deepen our understanding of e-h 



systems, leading toward the experimental observation. Com- 
pared with other systems undergoing condensation of bound 
pairs, electrons and holes in a quasi-thermal-equilibrium sys- 
tem in semiconductors have the following two notable charac- 
teristics: (i) They involve not only the e-h attractive Coulomb 
interaction but also the repulsive one between like particles 
(besides the Pauli exclusion principle), and (ii) they generally 
have different masses and mass anisotropics. To understand 
the coherent condensation of e-h pairs, therefore, we high- 
light these important characteristics in this work. Early and 
quite recent works^iiSSS based only on the BCS-like mean- 
field theory have shown that the mass difference and the mass 
anisotropy suppress the e-h BCS order. In the presence of such 
asymmetry, however, it is still unclear how the crossover from 
the e-h BCS state to the excitonic BEC evolves. In addition, 
little attention has been paid to the roles of electron-electron 
(e-e) and hole-hole (h-h) repulsive interactions in this prob- 
lem. 

In this paper, by calculating the transition temperature T c 
that directly reflects the excitonic BCS-BEC crossover, we 
clarify the effects of the repulsive interaction and the mass 
difference on the e-h pair condensation from the BCS to the 
BEC regime. A simple two-band Hubbard model with attrac- 
tive and repulsive on-site interactions is adopted to describe 
the e-h systems. Here we suppose that conduction electrons 
and valence holes, whose bands are isotropic, have infinite 
lifetime, and the number of electrons is equal to that of holes 
(N e = Nh). In our model the interaction strengths and the e-h 
density are treated as independent parameters. Several phys- 
ical quantities, such as the density of states, the density of 
occupied sites, and the quasiparticle weights, are calculated to 
discuss the properties of the normal phase. 

We employ the self-consistent f-matrix approximation 
(SCTMA) in our analysis. The SCTMA is an effective 
method for the BCS-BEC crossover problem at finite tem- 
peraturesr 4 ^ 5 ^ and is a conserving approximation based on 
the Baym-Kadanoff theory. This approximation deals cor- 
rectly with two-particle correlations. Thus it becomes asymp- 
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totically exact in the low-density limit. In addition to the 
SCTMA, we also use a local approximation (LA), which is 
justified in high spatial dimensions. The procedure is to ne- 
glect the momentum dependence of the self-energy and the 
vertex function. 26 It is in the same spirit as the k-averaged ap- 
proximation 27 or the dynamical mean- field theory (DMFT). 28 
As known in the DMFT literature, the LA itself (without other 
approximations) becomes exact in the limit of infinite spatial 
dimensions and a good approximation for three-dimensional 
systems. Although the SCTMA combined with the LA is not 
exact even in the infinite-coordination limit since only ladder 
diagrams are summed up, it has been successful for the su- 
perconductivity of the single-band attractive Hubbard model 
in high dimensions^ In particular, the successive interpola- 
tion between the BCS limit with T c oc exp(f /{/') and the BEC 
limit with T c °c t 2 /U' can be described well, 17,18 where U' 
and t denote the attractive interaction and the transfer energy, 
respectively. Therefore, we extend the scheme to our two- 
band model, and expect that our results are valid for three- 
dimensional bulk systems. 

This paper is organized as follows. In Sec. II, the SCTMA 
combined with the LA is applied to the normal phase for the 
e-h two-band Hubbard model. In Sec. Ill, several physical 
quantities are calculated and properties of the normal phase 
are discussed. In Sec. IV, the transition temperature from the 
normal phase to the e-h pair condensed phase is presented as a 
function of the e-h attraction strength for various parameters. 
We examine in Sec. IV A the effect of the e-e (h-h) repulsion 
on the transition temperature. In Sec. IV B, the effect of the 
mass difference on the transition temperature is analyzed and 
features of the BCS-BEC crossover in the e-h system are dis- 
cussed. Concluding remarks are given in Sec. V. 



II. FORMULATION 

We consider an e-h system described by the two-band Hub- 
bard model. The Hamiltonian is given by 



(ij).oa=e,h ja,a 
j,a jaa 1 



(1) 



where c"^ (c 1 ^) denotes a creation operator of an electron (a 
hole) with spin a = {T, |} at the ;'th site and n" a = c"^c" CT 
with a = {e,h}. The quantities t e (f/,) and \i e (jj,/,) are the 
transfer integral of the electrons (holes) between the nearest- 
neighbor sites and the chemical potential measured from the 
center of the bare electron (hole) band, respectively. The on- 
site Coulomb interaction of the e-e (h-h) repulsion and that of 
the e-h attraction are expressed by U and —U', respectively. 

We apply the SCTMA to the model Q. Feynman diagrams 
contributing to the self-energy of electrons and holes in the 
normal phase within the SCTMA are shown in Fig. 1, where 
all particle-hole (particle-particle) ladder diagrams are taken 
into account with respect to the interaction U (—!/'). The ex- 
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FIG. 1: Feynman diagrams for the self-energy in the normal phase. 
The solid line denotes the electron or hole Green's function G a 
where a = h(e) for a = e (h). The spin weight is 2 for the diagram 
of the vertex function r a a- 



plicit expression of the self-energy is given by 

Z a (k,(On) = ^E e ' (v "' +ffl " )0+ r aa (q,v m ) 
" q.v„, 

x G a (q + k,v,„ + ffl„) 

+ ^-E e, ' (v '"^" )0+r ««( < j' v ») 

N q-v,„ 
xGtt(q-k,v ffl -(i)„), 



(2) 



where TV is the total number of lattice sites, T denotes the tem- 
perature, co,, = (2n + l)7tT, and v„, = 2nmT with integer n 
and m. The symbol + in the convergence factor denotes a 
time infinitesimally later than T = 0. Here the spin index a 
is omitted because we now consider the spin-symmetric case, 
but the spin weight (a factor 2) is counted. The single-particle 
Matsubara Green's function G a (k, CO,,) defined by the Fourier 
transform of — (T-ccJ a (r)cJ, g ) is expressed in terms of the self- 
energy (|2ji as 



G a (k,0^) = 



icdn+Ha-eg -Z„(k,cow)' 



(3) 



where e£ is the band dispersion of the noninteracting elec- 
trons/holes. The two-particle vertex functions T aa i (q, v m ) are 
obtained as the infinite geometric series, 



r«a(q,V m ) = U + UK aa (q,V m )U- 
U 



(4) 



(5) 



l-UK a a(q,V m )' 
r a a(q,V m ) = -U'+(-U')K a& (q,V m )(-U')-\ 

-U' 

l + U'K aei (q,V m )' 
where the pair propagators are 
T 

K aa (q,v m ) = -— £ G a (k,d) n )G a (q+k,v m + a) n ), (6) 



k,<»„ 



K aa (q,V m ) = -— £ G a (k,ffl«)G a (q-k,V m -cO. (7) 
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By using the self-consistent set of Eqs. the 
model Q can be solved within the SCTMA. However, a 
straightforward calculation is very difficult especially for 
three-dimensional systems with large system size N. In this 
work, then, we combine a LA with the SCTMA to simplify 
the above self-consistent calculationi^kSI The procedure is 
performed by neglecting momentum dependence of the self- 
energy and the vertex function [i.e., E a (k, CO,,) — > T. a (co„) and 
r«a' (<L v »i) ~~ ¥ r aa / (v m )]. This analysis is equivalent to solv- 
ing an effective single-impurity problem of the DMFT within 
the SCTMA. By making the LA, the self-consistent Eqs. 10- 
can be written in terms of only local functions: 

Z a (<o n ) = T^e ,(v '" +m " )Q+ r aa (v m )G a (v m + co„) 

V m 

+ 2TY,e i(v '"- m » )0+ r a& (v m )G & (v m - toy,), (8) 
v m 

where the local Green's function is 



G a (co„) = -£G a (k,a)„) 



p2(«0 



iCOn + ^a-e-LaicOn)' 



(9) 



with the noninteracting density of states p„(e), and the local 
vertex functions are 



Laa(v m ) — 



u 



i-uK aa ( Vm y 
-u' 



(10) 
(11) 



l+U'Kaa{v m ) 
with the local pair propagators 

K aa (v m ) = -TY,G a (co n )G a (v m + co n ), (12) 

Kaa{v,„) = -TY,Ga(C0 n )Ga{v m -C0„). (13) 

m„ 

To determine the transition temperature of the conden- 
sation of e-h pairs, we examine the e-h singlet pair sus- 
ceptibility 29 x(q,V m ), defined by the Fourier transform of 
-(7V4_ CT (T)c5 (T (T)4 t CT (T')c5;_ (J (T')), which is given by the 
ladder terms corresponding to r e /,(q, V m ), i.e., 



z(q,v m ) = 



K e h(<l,V„,) 

l+U'K eh (q,v m )' 



(14) 



Here note that within the present SCTMA the repulsive in- 
teraction U does not appear explicitly in Eq. dl4> . but it in- 
fluences the e-h pair susceptibility through the single-particle 
Green's functions in Eq. (0. If the uniform static e-h pair 
susceptibility %(0,0) diverges for T \ T c , it is a signal of the 
onset of the e-h pair condensation (the Thouless criterion). 
Thus the transition temperature T c can be determined as the 
temperature satisfying the condition 



Since the self-energy is local (i.e., momentum independent), 
the pair propagator K e f,(0,0) in Eq. d!5t can be expressed in 
terms of only local functions as 



K eh (0,0) = --£G e (k,(H,)G fc (-k, 



-(On) 



k,co„ 



G e {co n )-yG h (co n )* 



(16) 



where Ca(co„) = ico„ + fl a - and y = t h /t e is de- 

fined as the e-h band mass ratio m e /m/,. Here we have as- 
sumed the relation e£ = ye£. In this paper we use p„(e) = 
-\/4? 2 — £ 2 /(27Tf 2 < ) as a typical density of states of three- 
dimensional systems. 



III. PROPERTIES OF NORMAL PHASE 

In this section, properties of the normal phase above T c are 
discussed by using the self-consistent solution of Eqs. 0- 
( II 3b . We analyze the interacting density of states, the density 
of sites occupied by electrons and holes, and the quasiparticle 
weight, for various values of £/, U', y — th/te, and the particle 
density n. In the following we consider the cases of n = 0.25 
and 0.5 (and n = 1 only in the next section), where the particle 
density is defined as 

n = ^{n%) = 2TY,e i<0 "° + G a (co n ), (17) 

a a>„ 

under the condition of charge neutrality (n e a ) = («^). 
Throughout the paper the quantity t e + 1„ is taken as the en- 
ergy unit. 



A. Density of states 



Figure[2]shows the interacting density of states per spin, 



Pa{(o) = lmG a (co + iO + ), 



(18) 



i+t/Xft(o,o)=o. 



(15) 



for various values of U' . Parameters are chosen as n = 0.25 
(1/8 filling), T = 0.04 (that is, above T c ), and y = 1 with the 
fixed U — (a) and (b) 2. Note that p e ((o) = Pi,{®) for y= 1. 
Here we have obtained the retarded Green's function G a {(0 + 
i0 + ) from the numerical analytic continuation by using the 
Pade approximation. 

When [7 = 0, p a ((o) exhibits a characteristic two-peak 
structure as U' is increased: a sharp peak (i.e., the quasi- 
particle coherent peak) appears at co ~ and its weight de- 
creases, while a broad incoherent peak develops at co ~ U'. 
We should remark that in the low-density limit (n — ► 0) an e-h 
bound state is formed for U' > 1 and the binding energy is 
given by £/, ~ U'. Hence, such a behavior of the density of 
states implies that local e-h pairs (excitons) tend to be formed, 
which will become clear also from the analysis of the density 
of occupied sites (Fig.|3|l- 
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FIG. 2: The density of states p«((o) for n = 0.25, T = 0.04, and 
7 = tf,/t e — 1 . All the energies are scaled by t e + f/, . 

When U = 2, p«(fi)) is broadened by U, so it is already 
strongly renormalized even for U' = 0. This indicates that U 
gives the dominant contribution to the band gap renormaliza- 
tion. For U' < U (U' = 1), in contrast to the case of U = 0, 
the broad incoherent peak is not observed. This means that the 
excitonic correlation is suppressed by U. For U' > £/ (U' = 2 
and 3), p a (co) shows the two-peak structure again. Compared 
with the case of U = 0, however, the coherent (incoherent) 
peak is enhanced (suppressed) especially for U ~ U', This 
is due to the competition between U and U', Similar behav- 
ior has been found more prominently in the DMFT for n = 1 
(half filling) There, the metallic phase is stabilized for 
U ~ U' < 5 between the Mott-Hubbard and biexcitonlike in- 
sulating phases. Except in such a case, the density of states in 
Fig-Elsuggests fundamentally that the metallic character tends 
to be lost because of both the weight shifting due to U' and the 
overall broadening due to U. 

Here we note that within the conserving approximation 
such as the present SCTMA, the system never becomes an in- 
sulator [pa(0) = 0] with no broken symmetry for any values 
of U and U 1 . This is a well-known failure of this scheme. 26 ' 32 



B. Density of occupied sites 

Examining the density of occupied sites is an effective way 
to see how the excitonic correlation increases in the normal 
phase. Now we focus on the density of sites quadruply oc- 
cupied by electrons and holes, = Y^aa'^ja^jc')- This * s 
local susceptibility at equal times and can be written in terms 
of the local pair propagator (tOl : 

We first note the limiting cases. For the noninteracting case 
D e h = n 2 - If a ll electrons and holes are perfectly bound as lo- 
cal e-h pairs by U', the density of occupied sites is expected 
to become D e f, = n for U > U' and = 2n for U < U', The 
former case corresponds to an exciton phase, and the latter 
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FIG. 3: The density of sites quadruply occupied by electrons and 
holes as a function of U 1 for n = 0.25, T = 0.04, and 7 = 1 (these 
parameters are the same as in Fig.|2j with fixed U = 0, 1, and 2. 



to a biexciton phase where two electrons and two holes (i.e., 
two excitons) gather on a site. Actually, exciton and biexciton 
insulating phases with the aforementioned features have been 
obtained by using DMFT combined with the exact diagonal- 
ization method, ^-i^ 

In Fig.|3j the U' dependence of D e /, is shown for n = 0.25, 
T = 0.04, and 7=1 (corresponding to Fig. |2j with fixed 
U = 0, 1, and 2. For U' = 0, the density of occupied sites 
becomes D e j t ~ n 2 = 0.0625, that is the value in the nonin- 
teracting case, regardless the value of U. With increasing U', 
D e h increases monotonically. This tendency clearly indicates 
the development of the excitonic correlation toward the for- 
mation of local e-h pairs. To detect trends more evidently, we 
drew the thin dotted line (D^ = n = 0.25) as a loose guide 
for a criterion of local e-h pair formation. Here, of course, we 
should keep in mind that the system is as metallic as ever, as 
seen in Fig. [2] even if is above the line. One can find that 
the excitonic correlation is suppressed by U for U ~U'. In 
the large-f/' region where the effect of U is suppressed, 
increases beyond 0.25 as a result of the enhancement of the 
biexcitonic correlation. However, D e /, does not reach 2n = 0.5 
corresponding to the biexciton phase. This is because the t- 
matrix approximation cannot deal adequately with the four- 
body correlation, irrespective of the LA. Hence problems of 
biexcitons and their condensation are out of the scope of this 
paper. 



C. Renormalization factor 

Here we consider the case that the electron and hole band 
masses are different (77^ 1). Evaluating the renormalization 
factor, we examine the properties of the normal phase in the 
case with mass difference. We define the renormalization fac- 
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FIG. 4: The renormalization factor Z a as a function of U for n = 0.5, 
T = 0.05, and y = 0.5 with fixed U' = 0,l, and 2. The solid (dashed) 
curve denotes Z e (Z/,). 



tor Z a by 



a!£ a (<»„) 



ia> n — >0 



1 



Im[E a (/7rr)-E a (~^r)] 
2nT 



(20) 



The quantity Z a at T — > has generally several physical mean- 
ings characterizing the Fermi liquid phase, e.g., the weight of 
the quasiparticle coherent peak in the density of states and the 
jump at the Fermi wave number in the momentum distribution 
function. Since the self-energy is momentum independent, it 
also gives directly the effective mass enhancement of quasi- 
particles, m a /m* a . 

In Fig. |4j we show the U dependence of Z a for n = 0.5 
(quarter filling), T = 0.05, and y = 0.5 (the hole band mass 
is twice that of the electron) with fixed U' = 0, 1, and 2. For 
U' = 0, both Z e and Z/, are monotonically decreasing func- 
tions with respect to U. Because of the narrow hole band, 
the effective repulsive interaction between holes U /(2th) is 
quite strong, leading to a large reduction of Z/, (dashed curve). 
The presence of U' (U' = 1 and 2) reduces both Z e and Z/,. 
However, the behaviors of Z e and Z/, for U' ^ are qualita- 
tively different from those for U' = 0: Z/, decreases slowly 
and monotonically with increasing U, while Z e has a slight 
hump (solid curve). The origin of this hump is the competi- 
tion between the effective repulsion U / (2t e ) and the effective 
e-h attraction U' /(t e + ?/,), which raises the mobility of elec- 
trons. When y = 0.5, therefore, the hump of Z e appears at 
U ~ 2t e /(t e + t h ) xU' = 2U'/(1 +y) ~ 1.3 and 2.7 for U' = 1 
and 2, respectively. For holes the influence of the competi- 
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FIG. 5: The transition temperature T c as a function of U' for n = 0.25 
and y = 1 with fixed U = 0, 1, and 2. The dotted line denotes the 
result from BCS theory. 



tion between U /(2f/,) and U' j[t e + ?/,) is quite weak, and Z/, 
is almost constant in the region of U < 2yU' / (1 + 7) — 0.7 
and 1.3 for U' = 1 and 2, respectively. Here we point out that 
Z a becomes insensitive to U for U' = 2. When U' = 2, the 
tendency to form bound e-h pairs begins to grow strongly, in- 
dicating very weak U dependence. This fact will also appear 
in the behavior of the transition temperature discussed in the 
next section. 



IV. TRANSITION TEMPERATURE 

In this section, based on the self-consistent solution and the 
condition defined in (1 1 51 . we evaluate the transition tempera- 
ture T c from the normal phase to the e-h pair condensed phase. 
By analyzing the U' dependence of T c for various U and y, the 
effects of the repulsive interaction and the mass difference on 
the transition temperature are discussed. 



A. Effect of U 

First, we consider the effect of the repulsive interaction U 
on the transition temperature T c . The calculations in this sub- 
section are restricted to the mass ratio y = 1 . 

Figure |3 shows the U' dependence of T c for n = 0.25 and 
y = 1 with fixed U = 0, 1, and 2. The e-h pair susceptibility in 
the normal phase diverges for T \ T c , which means that the 
e-h pair condensed phase is realized for T <T C . For U = 0, T c 
can be described well by the BCS result (dotted curve) in the 
weak-coupling region (U 1 < 0.3), that is, T c exp(— A/U') 
with a constant A. With increasing U' , the solid curve de- 
viates from the BCS result, reaches a maximum at U' ~ 1, 
and then decreases as l/U' for large U' . The U' dependence 
of T c for large U' is related to the behavior of the BEC tem- 
perature of a lattice boson system 17,18 - 26 with kinetic and po- 
tential energies of order 1 /U 1 ? A Obviously this result implies 
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FIG. 6: The transition temperature T c as a function of U' for y = 1 
and U = 2 with fixed n = 0.25, 0.5, and 1.0. The inset shows the case 
of U = for n = 0.25 (solid line), 0.5 (dotted line), and 1.0 (dashed 
line). 



the BCS-BEC crossover, as expected. Note that the ladder 
terms for the dressed Green's function is essential to extract 
the above successive crossover. 24 In the presence of U (U = 1 
and 2), one can see the reduction of T c for U' < U . This fact 
can be understood as a consequence of the suppression of the 
excitonic correlation discussed in Figs. [2] and [5] In contrast, 
T c is largely insensitive to the repulsion U for U' > U . This 
supports the validity of physical picture of strongly bound e-h 
pairs that behave almost like neutral bosons for U' > U, which 
is consistent with the analysis of the density of occupied sites 
in Fig. |3l 

Next, the effect of the repulsion U on the transition tem- 
perature is examined by changing the e-h particle density. In 
Fig. |6j the U' dependence of T c with y = 1 and U = 2 is 
shown for n = 0.25, 0.5, and 1 (half filling). The inset dis- 
plays the corresponding U' dependence of T c for U = 0. When 
U = (inset of Fig. [6}, the transition temperature T c increases 
slightly in the weak-coupling region (U' <C 1) as the particle 
density n is increased. This behavior is quite reasonable for 
U' <C 1 from the following reason. Within the BCS theory the 
effective attractive interaction between electrons and holes is 
roughly given by p^(e")t/' where e" is the Fermi energy for 
an uncorrected system [see Eq. (BU I. Hence the effective at- 
traction reaches maximum when n = 1, i.e., the band is half 
filled, in the case of the semicircular density of states. The 
presence of U, however, completely changes this tendency. 
As seen in the main figure of Fig.[6]([/ = 2), T c is sufficiently 
suppressed by U for U' < U as n approaches 1. When we 
simply consider the effect of U in the BCS regime, the ef- 
fective interaction p^(e")C/' could be replaced by p a (Sp)U', 
where Pa(e") is strongly renormalized by U , as seen in Fig. [2] 
The renormalization becomes strong as n approaches 1 M> As 
a result, the effective interaction is reduced by U as n — ► 1, 
leading to the suppression of T c . Meanwhile, for large U' the 
contribution of U to T c is very small regardless of the value of 
«. 
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FIG. 7: The transition temperature T c as a function of U' for n = 0.25 
and U = U' with fixed 7=1, 0.5, and 0.2. In the inset, T c is shown 
as a function of U'/(t e th) for 7 = 1 (solid line), 0.5 (dotted line), 
and 0.2 (dashed line), where the decimal logarithm is taken for both 
axes. The dot-dashed line represents the exact result 1221 in the limit 
of U — > 00 and large U' . 

B. Effect of mass difference 

Finally, we investigate the effect of the mass difference 7 
on the transition temperature T c . 

We plot the U' dependence of T c for n = 0.25 and U = U' 
in Fig.0 with fixed 7=1, 0.5, and 0.2. The mass difference 
has the effect of reducing the transition temperature, which is 
qualitatively consistent with the result of the BCS-like pairing 
theory by Mizoo et aljoL In addition to the reduction of T c in 
the weak-coupling BCS regime, our result also involves the 
suppression of T c in the strong-coupling BEC regime. Such 
a behavior of T c can be understood analytically in the weak- 
and strong-coupling limits as follows. 

In the weak-coupling limit, the transition temperature can 
be estimated by using the Green's function including only the 
Hartree term in its own self-energy, which is given by 

*=-i.i3^-K-s^). (2l) 

where w" is a cutoff energy of order e". Note that t e p ( J(ep ) — 
ti t p® {e'p). Equation i2H indicates that the effective e-h attrac- 
tive interaction is given by f a p^(£")J/', as mentioned above. 
Moreover, another important point derived from Eq. (12 It is 
that the system is characterized by an energy scale t e + th de- 
termining the dimensionless effective coupling strength in the 
weak-coupling region. In our model the quantity (r e + ?/?) ' is 
proportional to m e m^/(m e + namely, the reduced mass. 
Therefore, the transition to the e-h pair condensed phase in 
the weak-coupling region is related to the relative motion be- 
tween electrons and holes, implying the BCS regime. The 7 
dependence of T c in the BCS regime can be roughly evaluated 
from that of the cutoff energy w", where we should recall that 
the argument of the exponential in Eq. (12 It does not depend 
on the mass ratio 7 when t e + f/, is taken as the energy unit. 
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The cutoff energy of holes should be yw e c (~ yt e ) since 
Bp = Y £ f (~ y { e)- Thus the coefficient of Eq. d2Tl becomes 
a/ w£w* ~ , and so it is found that T c in the BCS regime 
would be approximately proportional to -Jy/ (1 + y). 

On the other hand, in the strong-coupling limit, the 
model Q can be mapped onto a single-band attractive Hub- 
bard model with the interaction —U' and the hopping t a at 
JJ = oo. For large U' , this attractive Hubbard model can map to 
a hard core Boson model with the kinetic energy 2t e tf,/U' and 
the potential energy (t\ +tf l )/U'. Using the standard mean- 
field theory, which becomes exact in the limit of the large 
coordination number of the lattice, we can obtain the BEC 
temperature in the limit of U — > °° and large U' as 35 

t^bec _ 2n - 1 

U' ln[n/(l-n)]' k ' 

This expression indicates that the system is characterized 
by an energy scale t e tf,/(t e + f/ 7 ) determining the dimension- 
less effective coupling strength in the strong-coupling region, 
which is related to the motion of the center of mass since 
(t e + th) I (t e th) m e +nif l . The U' dependence of T c in the 
strong-coupling region of Fig.[7]can be described asymptot- 
ically in terms of t e tf,/U'. In the inset of Fig.Q the transi- 
tion temperature is shown together with the result of Eq. 122k 
where U'/(t e tf,) is taken as the x axis and the decimal loga- 
rithm is taken for both axes. The transition temperature T c for 
various mass ratio y tends to behave linearly with the slope of 
— 1 (i.e., T c oc 1 /{/') and merges in the large-f/' region, imply- 
ing the BEC regime. As seen from Eq. (I22L the / dependence 
of T c in the BEC regime is given by T c « t e tf, = y/(l + y) 2 . 
Therefore, it is confirmed that for the fixed U' and n the tran- 
sition temperature is suppressed by the e-h mass difference in 
the BEC regime, as well as in the BCS regime. 

Viewed in this way, the characteristic energy scale and the 
e-h mass ratio dependence of the transition temperature be- 
come clear in both the BCS and BEC regimes. We especially 
emphasize the point that the excitonic BCS-BEC crossover 
can be marked by the change of the characteristic energy scale 
from t e +t/ r to t e th/(t e +th), which is not necessarily evident 
in the standard scenario 19 of the BCS-BEC crossover. 



V. CONCLUDING REMARKS 

In this paper, by applying the self-consistent f-matrix ap- 
proximation with the LA to the e-h two-band Hubbard model, 
we have discussed the properties of the normal phase and the 
transition temperature to the e-h pair condensed phase. In the 
analysis, the effects of the repulsive interaction between like 
particles and the mass e-h difference have been given special 
attention. 

From the behavior of the calculated physical quantities, 
we have found that the development of the excitonic correla- 
tion induced by the e-h attractive interaction is suppressed in 
the presence of the repulsion between like particles and also 
the mass difference. In particular, it is noteworthy that the 
competition between the effective interactions U / (2t a ) and 
U'/(t e + f/j) play the key roles in controlling the formation of 



excitonic bound pairs in the normal phase. This is a remark- 
able feature of the two-band e-h system. 

The behavior of the transition temperature has shown that 
the present SCTMA combined with the LA can describe well 
the excitonic BCS-BEC crossover at finite temperatures. In 
the BCS regime of weak e-h attraction, the repulsive inter- 
action sufficiently suppresses the transition temperature for 
U > U'. In contrast, it does not affect the transition temper- 
ature for U < U' in the BEC regime of strong e-h attraction. 
These results imply that the competition between the effec- 
tive repulsion and attraction plays an important role in under- 
standing the properties of the e-h pair condensed phase and 
the crossover behavior also below T c . Furthermore, we have 
found that in the whole BCS-BEC regime the e-h mass differ- 
ence leads to large suppression of the transition temperature. 
This analysis in the case of y ^ 1 has allowed us to capture the 
BCS-BEC crossover as the change of the characteristic energy 
scale from t e +tf, to t e tf,/(t e + f/j), where the former is related 
to the relative motion and the latter to the motion of the center 
of mass. 

Let us now discuss limitations of the present SCTMA. The 
n dependence of T c in the large-f/' region is opposite to that 
of the exact result in the limit of U — > °° and large U' given by 
Eq. J22I . The result of Fig. [6] suggests an decreasing function 
of « while the exact result is an increasing function of n. It is 
not clear at present why these are not consistent even qualita- 
tively. At least we should remind ourselves that the f-matrix 
approximation is valid for the low-density limit. 

Another difficulty is that the SCTMA cannot describe in- 
sulating states without any symmetry breakingr^ 3 ^ such as 
the exciton and biexciton phases. Here we should remark that 
the SCTMA for the single-band attractive Hubbard model 26 
does not yield the pseudogap behavior (or the pairing transi- 
tion 3 ^ 3 ^ 3 ^ 3 ^) as seen in high-r f superconductors. The pseu- 
dogap phase in the context of superconductors is deeply re- 
lated to the exciton phase in e-h systems where incoherent e-h 
bound pairs (which do not condense but are insulating) are 
formed. In our model, the DMFT combined with the exact 
diagonalization method^L 3 ^ indicates the appearance of exci- 
ton and biexciton insulating phases that already show up for 
U' ~ 2. Nevertheless, it is still notable that even within the 
SCTMA one can see the tendency to form e-h bound pairs 
above T c by calculating the density of states and the density of 
occupied sites. However, a unified theory for the exciton Mott 
transition and the excitonic condensation would requi re non- 
perturbative approaches such as the full DMFT (Refs. l3lj33l 
and|39) and Monte Carlo methods.!^ 

The long-range part of the Coulomb interaction was not 
considered here. It may become crucial particularly in the 
BEC regime, which motivates future work. 

As has been argued in Refs. ElandE! clarifying the BCS- 
BEC crossover problem from the optical response attracts 
great interest and is important from not only the theoretical 
but also the experimental aspect. For this purpose, the present 
SCTMA will be extended to the e-h pair condensed phase be- 
low T c . The analysis for the properties of the condensed phase 
and optical response will be reported soon in a forthcoming 
presentation^ 
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